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Hughes  are  also  included. 
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™  1^  SHALLOW-WATER  SPECTRUM 
DESCRIPTION  AND  APPLICATIONS 

PART  I:  INTRODUCTION 

1.  Determination  of  realistic  estimates  for  wave  conditions  in  shallow 
water  has  always  been  important  in  coastal  engineering.  Initially  monochro¬ 
matic  wave  theory  was  used  exclusively,  but  it  is  ..ow  generally  accepted  that 
an  irregular  wave  approach  is  preferred.  The  purposes  of  this  report  are  to 
summarize  recent  advances  in  the  specification  of  a  shallow-water  self-similar 
spectral  form,  to  discuss  the  range  of  practical  usage,  and  to  provide  engi¬ 
neering  examples  using  the  new  formulations. 

2.  The  first  quantitative  predictions  of  ocean  waves  evolved  from  ob¬ 
servations  of  waves  at  a  point  in  space  and  correlations  between  such  parame¬ 
ters  as  mean  wave  height,  significant  wave  height,  and  mean  wave  period.  In 
the  mid-1950's  the  spectral  description  of  the  sea  state  was  forwarded.  This 
description  allowed  for  more  detail  by  distinguishing  between  the  different 
frequency  bands.  Phillips  (1958)  suggested  that  there  should  be  a  region  of 
the  spectrum  of  wind-generated  deepwater  gravity  waves  in  which  the  wave 
energy  density  has  an  upper  bound  given  by  the  following  expression: 

E  (f)  =  ag^f"^(2n)"^  (1) 

m 

where  f  is  frequency,*  g  is  gravitational  acceleration,  and  a  was  as- 

-3 

sumed  to  be  a  universal  constant,  approximately  8  x  lO  .  This  region  of  the 
spectrum  is  called  the  equilibrium  range.  In  scalar  wave  number  space  this 
equilibrium  range  becomes 


F(k)  =  (2) 

where  p  is  a  constant  and  k  is  the  wave  number  equal  to  2n/L  when  L  is 
the  wave  length.  The  limit,  defined  by  Equation  1  or  2,  was  thought  to  be  a 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  A) . 
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result  of  a  limiting  wave  steepness  at  each  frequency  beyond  which  deepwater 
wave  breaking  would  occur.  In  other  words,  any  additional  energy  input  in  the 
spectrum  at  that  frequency  would  result  in  waves  breaking  and  a  transfer  of 
wave  energy  from  that  frequency  through  the  mechanisms  of  dissipation  and 
wave-wave  interactions.  This  represents  the  equilibrium  or  steady-state  wind 
sea  situation  in  deep  water,  and  Equation  1  describes  the  portion  of  the  spec¬ 
trum  to  the  high  frequency  side  of  the  single  spectral  peak. 

3.  A  spectral  shape  for  fully  developed  deepwater  waves  was  advanced  by 
Pierson  and  Moskowitz  (1964)  containing  Phillips'  equilibrium  range 
(Equation  1).  It  is  expressed  as 

-5/4(f/f  )'^ 

E  (f)  =  E  (f)e  (3) 

pm  m 

where  f^  is  the  frequency  of  the  spectral  peak.  The  additional  exponential 
term  provides  the  low  frequency  forward  face  of  the  spectrum  and  a  broad, 
smooth  peak  region.  Equation  3  is  based  on  extensive  field  data.  Estimation 
of  f^  for  fully  developed  seas  (i.e.,  no  further  wave  growth  could  occur  at 
that  windspeed)  was  empirically  determined  in  terms  of  10-m-high  windspeed  U 
as 


f  =  0^ 

m  2nU 


(4) 


Integration  of  Equation  3  provides  the  total  energy  E  in  the  spectrum. 

1/2 

Longuet-Higgins  (1952)  has  shown  that  4  times  the  variance,  or  4(E)  ,  pro¬ 

vides  a  close  estimate  of  significant  wave  height  in  deep  water.  Thus  a 
method  exists  for  hindcasting  and  forecasting  fully  developed  deepwater  wave 
heights  and  periods. 

4.  One  drawback  to  the  fully  developed  representation  was  that  at 
higher  windspeeds  the  wind  seldom  held  steady  for  the  length  of  time  necessary 
for  fully  saturated  seas  to  develop.  In  other  cases  the  fetch  over  which  the 
wind  acted  wasn't  long  enough  for  fully  developed  conditions. 

5.  The  data  from  the  Joint  North  Sea  Wave  Program  (JONSWAP)  were  used 
by  Hasselmann,  et  al.  (1973)  to  extend  the  Pierson-Moskowitz  spectral  density 
representation  (Equation  3)  to  partially  developed  wave  conditions  by  the  ad¬ 
dition  of  another  factor: 
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where 


“  =  0  076 


new  factor 

-0.22 


-0.33 


■•(S 


-0.143 


(Mitsuyasu  1981) 


U  =  windspeed 


X  = 

fetch  distance 

(a 

=  0.07  for 

f 

>  f 

m 

K 

=0.09  for 

f 

<  f 

\  b 

m 

6. 

The  effect 

of  ' 

the  1 

spectra  which  are  typical  of  growing  wind  seas  in  deep  water.  Correlations 
were  used  to  arrive  at  the  parametric  relationships  given  by  Equations  6-8. 
Other  self-similar  spectral  representations  have  been  proposed  for  deep  water 
(Toba  1973  and  Kruseman  1976),  but  the  JONSWAP  form  is  the  most  widely  known. 


PART  II:  SHALLOW- WATER  DEVELOPMENTS 


Equilibrium  Range 

7.  Kitaigorodskii  et  al.  (1975)  examined  the  possibility  that  an  equi¬ 
librium  range  existed  also  for  the  finite  depth  horizontal  bottom  case.  They 
noted  that  while  both  forms  of  the  equilibrium  range  given  for  deep  water 
(Equations  1  and  2)  are  valid,  only  one  of  them  could  hold  simultaneously  for 
deep  and  shallow  water.  Since  wave  breaking  occurs  in  both  deep  and  shallow 
water,  it  is  reasonable  to  believe  that  the  similarity  form  must  be  consistent 
between  deep  and  shallow  water.  Field  observations  led  Kitaigorodskii  et  al. 
to  conclude  that  the  wave  number  expression  (Equation  2)  was  the  valid  scaling 
for  both  the  deepwater  and  the  shallow-water  wind  wave  equilibrium  ranges. 

This  has  been  confirmed  by  other  investigators  (Gadzhiyev  and  Kratsitsky  1978, 
Vincent  1982). 

8.  Taking  the  shallow-water  limit  of  the  linear  wave  dispersion  rela¬ 
tion,  i.e.  , 

lim  r(27rf)^  =  gk  tanh  (kh)l  or  2nf  =  (gh)^^^k  (9) 

h^o  L  J 


_3 

where  h  is  the  water  depth,  it  is  seen  that  a  k  •  dependence  is  equivalent 
-3 

to  an  f  dependence  in  shallow  water.  Thus,  in  frequency  notation,  the 

slope  of  the  high  frequency  side  of  the  nearly  saturated  wind  wave  spectrum 

-5  -3 

transforms  from  an  f  to  an  f  during  the  shoaling  process  for  the  part 

of  the  spectrum  that  can  be  described  by  linear  wave  theory. 

9.  Kitaigorodskii  et  al.  (1975)  developed  a  frequency  dependent  factor 

<l>(27tf,h)  to  transform  the  f  ^  deepwater  range  factor  E  (f)  to  its  finite 

m 

depth  water  equivalent,  resulting  in  the  following  expression  for  the  finite 
water  depth  equilibrium  range: 


where 


and  w  =  2nf. 


E^(f,h)  =  E^(f)  •  «l>(27tf,h) 

-  ag^f"^(2n)'^  •  <J>(27tf,h) 


<t>(27Tf  ,h) 


k  (u»,h) 


k  ^(u),<») 


9k(iu.h) 

9u) 

9k(u).w) 

9u) 


(10) 


(11) 
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An  iterative  expression  arises  for  <t>(2nf,h)  when  linear  wave  theory  is  used 
in  Equation  11.  One  form  is  given  by 


with 


«l>(2nf,h) 


1-1 


sinh 


(12) 


(13) 


and  obtained  from  the  iterative  solution  of 

1 


R(u).  )  tanh 
h 


(14) 


The  function  4>  approaches  the  value  of  one  in  deep  water  and  a  value  of 
zero  as  thr  depth  decreases,  as  shown  in  Figure  1.  A  simple  approximation 
for  <l>(2nf,h)  is  given  in  Thompson  and  Vincent  (1983)  as 


<l>(2nf,h) 


U),.  <  1 


U).  >  1 

h 


(15) 


This  approximation  deviates  from  Equation  12  by  a  maximum  of  4  percent  in  the 

vicinity  of  u)j^  =  1.0  with  the  error  less  than  1  percent  outside  the  range 

0.8  <  u),,  <  1.3  . 
h 


Figure  1.  Kitaigorodskii  et  al.’s  «!>  as  a  function  of  U)j^ 

(Bouws  et  al.  1985a) 
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10.  Vincent  (1982)  pointed  out  that  the  estimate  of  the  upper  limit  on 
energy  density  in  shallow  water,  as  provided  by  Equation  10,  could  be  inte¬ 
grated  to  get  total  energy  if  a  low-frequency  cutoff  value  is  known,  i.e., 

00 

E  =  y*  ag^f‘^(27t)"^«J>(2nf,h)  df  (16) 

f 

c 


For  shallow  water,  when  the  primary  frequency  components  fall  in  the  range  of 
lu^  <  1  ,  and  assuming  that  f^  =  0.9f^  ,  Equation  16  becomes 


J)  dt  =  f  f-3  df 


when  Equation  15  is  used  to  approximate  <t>(2nf,h)  .  Integrating  Equation  17 
and  applying  the  limits  yields 


4(271) 


Using  the  energy-based  definition  for  significant  wave  height. 


H  =  4(E)' 
mo 


Equation  18  provides  an  estimate  of  the  depth-limited  significant  wave  height 


^  (agh)^/^T 

7t  m 


where  T  =  1/f 
m  m 

11.  Initially,  a  in  Equation  20  was  taken  to  be  Phillips'  constant 

equal  to  0.0081,  but  more  recent  results  have  provided  an  expression  for  a 

which  varies  according  to  depth,  windspeed,  and  peak  frequency  f^  .  This 

will  be  covered  shortly  in  more  detail.  Engineering  use  of  Equation  20  should 

be  limited  to  water  depths  between  the  surf  zone  out  to  depths  where  the 

primary  energy-containing  frequencies  are  defined  by  u),  =  1  . 

h 


TMA  Spectrum 


Self-similar  spectral  form 

12.  Bouws  et  al.  (1985a)  hypothesized  that  a  first  approximation  for  a 

finite  depth  wind  sea  spectral  shape  would  be  obtained  simply  by  substitution 
of  Kitaigorodskii  et  al.*s  expression  for  the  equilibrium  range  (given  by 
Equation  10)  for  the  deepwater  equilibrium  range  the  JONSWAP  Equa¬ 

tion  (given  in  Equation  5).  They  obtained  the  following: 

o  -s  -U  -5/4(f/f  )‘^  exp  [-(f/f  -l)^/2a^l 

ETfjA(f»h)  =  ag'^f  ^(2n)  >(27if,h)e  ”  y  L  ”  J  (21) 

with  <l>(2nf,h)  given  by  Equation  11. 

13.  Bouws  et  al.  named  their  self-similar  finite  water  depth  spectral 

shape  the  lilA  spectrum  by  combining  the  first  three  letters  of  the  three  data 

sets  used  for  field  verification  (Texel,  MARSEN,  and  ARSLOE),  each  of  which  is 

explained  below.  Equation  21  has  the  four  JONSWAP  parameters,  a  ,  y  ,  f  , 

in 

and  a  plus  the  depth  h  .  An  example  of  the  finite  depth  effect  of  the  THA 
spectral  form  is  shown  in  Figure  2  in  which  all  the  parameters  are  held  con¬ 
stant  except  the  depth.  Figure  3  is  the  same  presentation  but  plotted  in  wave 

-3 

number  space,  which  shows  the  k  scaling  from  deep  water  to  shallow  water. 
Field  verification 

14.  Bouws  et  al.  (1985a)  used  field  data  from  three  separate  studies  on 
shallow-water  wind  wave  growth  to  investigate  the  parameters  used  in  the  TMA 
spectral  representation.  MARSEN  (the  Marine  Remote  Sensing  Experiment  at  the 
North  Sea)  and  ARSLOE  (the  Atlantic  Remote  Sensing  Land-Ocean  Experiment)  were 
both  comprehensive  experiments  in  which  ocean  wave  measurements  were  but  a 
part  of  the  entire  program.  The  experimental  sites  are  both  on  the  continental 
shelf  with  depths  up  to  40  m;  but  the  ARSLOE  site  was  open  to  the  Atlantic 
Ocean,  while  the  MARSEN  site  was  located  in  the  southern  half  of  the  North  Sea. 
The  Texel  data  set  is  comprised  of  a  series  of  measurements  made  near  the  Texel 
lightship  west  of  Rotterdam  during  a  longlasting  northwesterly  storm  in  the 
central  and  southern  North  Sea. 

15.  The  combined  data  represent  conditions  with  windspeeds  ranging 
between  4  and  25  m/sec,  bottom  materials  ranging  from  fine  to  coarse  sands, 
bottom  slopes  ranging  from  1:150  to  nearly  flat,  and  depths  from  about  5  m 
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Figure  2.  A  family  of  wind  wave 
spectra  with  identical  JONSWAP 
parameters  in  frequency  space 
(TMA  spectral  form) 


Figure  3.  A  family  of  wind  wave 
spectra  with  identical  JONSWAP 
parameters  in  wave  number  space 
(THA  spectral  form) 


to  45  m.  An  extensive  description  of  the  data  set  is  given  in  Bouws  et  al. 
(1985a). 

16.  Bouws  et  al.  (1985a)  fitted  the  TMA  spectral  form  to  over  2,800 
wind  sea  spectra  to  test  its  viability  and  to  determine  if  any  parametric  re¬ 
lationships  could  be  established  linking  the  spectral  parameters  to  the  ex¬ 
ternal  wind  field.  In  general  the  fit  of  the  spectrum  was  of  the  same  quality 
as  the  fit  of  the  JONSWAP  spectrum  to  deepwater  wind  sea  spectra.  Examples  of 
the  fit  to  field  spectra  are  given  in  Figure  4. 


TMA  parametric  relationships 


17.  The  relationships  between  the  TMA  spectral  parameters  and  various 
dimensionless  quantities  such  as  dimensionless  peak  frequency  and  wave  number 
were  examined  by  Bouws  et  al.  (1985a).  They  determined  that  a  and  y  could 
be  expressed  by  the  following  empirical  expressions  for  all  water  depths: 


0  =  0.0078k 


Y  =  2.47k 


fmmrn 


fS?; 


US' 


630 

80/10/25.  8.35 
h  =  18  m 


where 


with 


(24) 


U  =  windspeed  at  10  m  elevation 

g  =  acceleration  of  gravity 

k  =  2n/L  =  wave  number  for  waves  at  peak  frequency 
m  m 

L  =  wavelength  associated  with  peak  frequency  f  from  linear  wave 
theory 


Examination  of  the  variation  of  o  in  the  fitting  of  the  spectra  to  the  TMA 
representation  gave  no  statistical  correlation,  so  it  was  suggested  that  the 
mean  values  of 


a 


0.07 

0.09 


f  <  f 
-  m 


f  >  f 


m 


(25) 


are  sufficient  when  using  the  TMA  spectrum. 

18.  A  certain  amount  of  scatter  was  present  in  the  plots  from  which  the 
parametric  relationships  were  derived.  Bouws  et  al.  (1985b)  demonstrated  that 
the  scatter  could  be  simulated  by  assuming  that  the  error  in  determination  of 
windspeed  was  normally  distributed  with  a  standard  deviation  between  10  and 

20  percent.  They  concluded  that  the  methods  used  to  determine  windspeed  could 
easily  contain  errors  on  that  order.  The  most  remarkable  consequences  of  the 
TMA  research  were  (a)  showing  that  the  effect  of  depth  could  be  incorporated 
solely  through  the  dependence  of  K  on  depth,  and  (b)  determining  that  there 
was  no  apparent  systematic  dependence  on  bottom  type. 

19.  The  preceding  equations  provide  a  means  of  predicting  the  limiting 

wind  sea  spectral  shape  in  finite  depth  waters  in  terms  of  the  windspeed  and 
the  peak  frequency.  The  procedure  is  to  determine  K  from  Equation  24,  find 
a  and  y  (Equations  22  and  23,  respectively),  and  then  calculate  E^j^(f,h) 
values  using  Equation  21.  The  <l>  function  is  determined  for  each  frequency 
using  either  Equation  12  or  (if  a  slightly  less  accurate  answer  can  be  toler¬ 
ated)  Equation  15.  Significant  wave  height  can  be  determined  by  numerically 
integrating  the  resulting  shallow-water  spectrum  and  then  applying  the  rela¬ 
tionship  for  given  by  Equation  19,  but  an  easier  approach  is  to  use 

Vincent's  depth-limited  significant  wave  height  expression  given  by 


Equation  20,  provided  that  tOj^  <  1  for  most  of  the  energy  containing 
frequencies . 

TMA  internal  parameterizations 

20.  Starting  with  Kitaigorodskii  et  al.'s  (1975)  basic  equation  for  the 
equilibrium  range  (Equation  2),  Vincent  (1985a)  was  able  to  derive  an  expres¬ 
sion  for  a  in  terms  of  significant  wave  steepness  e  .  This  expression  is 


the  following: 


2  2 

a  =  l6n  e 


where 


1/2 

e  =  significant  wave  steepness  =  E  /L 


E  =  total  energy  in  the  spectrum 

=  wavelength  associated  with  f^  as  defined  in  linear  theory 
This  result  is  identical  with  that  obtained  by  Huang  et  al.  (1981)  for  a 
JONSWAP  spectrum  in  deep  water.  Field  data  support  Vincent's  formulation; 
thus,  Equation  26  holds  for  all  water  depths. 

21.  Equating  Equation  26  and  Equation  22  gives  an  expression  for  K  in 
terms  of  e  which  can  then  be  put  into  Equation  23  to  yield 


y  =  66l4e 


which  should  be  valid  for  wind  seas.  The  main  advantage  of  this  internal 

parameterization  is  the  ability  to  specify  the  equilibrium  wind  sea  spectrum 

associated  with  a  given  energy  level. 

22.  It  is  an  interesting  exercise  to  show  that  Vincent's  depth-limited 

significant  wave  height  equation  (Equation  20)  and  the  expression  for  a 

(Equation  26)  are  the  same.  The  energy  E  in  Equation  27  can  be  given  in 

terms  of  H  from  the  relationship  of  Equation  19.  Substitution  into  Equa- 
mo 

tion  26  yields 


H  =  i  (a)^^^L 
mo  71  m 


(30) 


The  only  difference  between  Equation  30  and  Equation  20  is  the  1.1  factor 
which  Vincent  (1982)  arrived  at  by  assuming  a  cutoff  frequency  of  0.9  f^  . 
Equation  29  is  a  more  general  form  of  Equation  30,  and  it  can  be  applied 
throughout  the  shallow-  and  intermediate-water  depth  ranges  when  the  appropri¬ 
ate  expression  for  L  is  used  and  a  is  found  from  Equation  22. 

m 

23.  To  reiterate,  significant  wave  heights  H  found  using  either 

mo 

Equation  29  or  Equation  30  are  the  maximum  values  that  can  occur  at  that  depth 
for  the  given  values  of  peak  period  and  windspeed.  In  addition,  calcu¬ 

lated  for  very  shallow  water  can  be  substantially  less  than  the  significant 
wave  height  defined  as  the  average  of  the  highest  1/3  waves.  pro¬ 

vides  an  estimate  of  the  energy  contained  in  the  wave  spectrum,  while  H  pro- 

s 

vides  a  statistical  estimation  of  the  wave  heights.  This  variation  in  shallow 

water  is  due  to  the  increasing  nonlinearity  in  the  waves.  A  full  treatment  of 

this  topic  is  provided  in  Thompson  and  Vincent  (1985),  and  provision  is  made 

for  this  effect  in  a  later  section  of  this  report. 

2A.  Since  Equation  29  reduces  to  Vincent's  depth-limited  result  in 

shallow  water,  it  might  be  suspected  that  a  similar  exercise  could  be  performed 

for  deep  water.  Replacing  L  in  Equation  29  with  the  deepwater  relationship, 
2 

L  =  (g/27T)T  ,  replacing  a  with  Equation  22,  and  nondimensionalizing  re- 

mo  m 


suits  in 


=  0.0112 


when  the  exponent  of  K  in  Equation  22  is  taken  as  1/2.  A  similar  result  is 
found  if  the  deepwater  fetch-limited  equations  for  wave  height  and  wave  period 
in  the  Shore  Protection  Manual  (SPM)  (1984)  are  combined,  the  only  difference 
being  a  coefficient  of  0.0105  instead  of  0.0112. 

25.  The  expression  for  the  fully  developed  wind  sea  significant  wave 
height  in  deep  water  can  now  be  derived  by  replacing  T^  in  Equation  31  with 
the  fully  developed  wave  period  as  defined  in  Equation  4.  This  results  in 

gH 

— ^  =  0.238  (32) 


The  only  difference  between  Equation  32  and  the  expression  for  fully  developed 

deep  water  H  given  in  the  SPM  (1984)  is  the  value  of  0.243  instead  of 
mo 

0.238.  This  is  remarkable  because  the  SPM  equation  was  arrived  at  by 


consideration  of  the  time  necessary  to  reach  full  development  of  the  JONSWAP 
equations . 

26.  This  exercise  has  illustrated  that  the  internal  and  external 
parameterizations  for  a  can  be  applied  for  all  water  depths  outside  of  the 
surf  zone.  The  ability  to  recover  deepwater  results  is  particularly  interest¬ 
ing  since  the  parameterization  of  a  as  a  function  of  f  and  windspeed  was 

m 

done  using  data  largely  from  shallow-  and  intermediate-water  depths. 

Wave  Growth  Limitation 


Wave  period  limit 

27.  Vincent  and  Hughes  (1985)  examined  the  question  of  whether  or  not 
finite  depth  wind  waves  reach  a  fully  developed  state  where  the  peak  spectral 
frequency  becomes  fixed  and  no  longer  migrates  toward  lower  frequencies, 
regardless  of  the  energy  input.  Evidence  that  such  a  cutoff  frequency  may 
exist  comes  from  the  Lake  Okeechobee  and  the  Gulf  of  Mexico  data  used  by 


Bretschneider  (1958)  to  develop  shallow-water  wave  growth  curves.  Bret- 


schneider  plotted  versus  gh/U  and  determined  an  empirical  relation¬ 


ship  for  the  cutoff  frequency  in  shallow  water.  Vincent  and  Hughes  pointed 
out  that  the  quenching  of  atmospheric  input  of  energy,  which  occurs  in  deep 
water  when  the  celerity  of  the  waves  exceed  the  windspeed,  did  not  appear  to 
occur  in  shallow  water.  This  is  because  at  a  given  depth  all  waves  of  fre¬ 
quency  less  than  f*  have  the  same  celerity,  if  f*  is  the  highest  frequency 
where  the  shallow-water  dispersion  relation  is  a  valid  approximation.  If  the 
deepwater  growth  mechanism  is  assumed,  it  would  be  possible  to  have  windspeeds 
in  excess  of  the  wave  speed  and  still  have  wave  growth. 

28.  Vincent  and  Hughes  (1985)  proposed  that  the  cutoff  frequency  where 
shallow-water  wave  growth  would  stop  should  be  determined  by  a  frequency  where 
the  shallow-water  dispersion  relation  held  plus  a  further  shift  to  slightly 
lower  frequencies  due  to  resonant  wave-wave  interactions.  This  frequency  is 
defined  as 


u).  =  2nf 

hm  m 


(i) 


1/2 


=  0.9 


(33) 


and  it  arises  from  Kitaigorodskii  et  al.'s  (1975)  observation  that  the  tran- 
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sition  from  a  deepwater  f  equilibrium  range  to  an  f  range  for 
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shallow-water  waves  occurred  when  u)j^  :s;  1  .  The  slight  downward  shift  in 
frequency  is  empirically  represented  by  the  0.9  value. 

29.  Rearranging  Equation  33,  multiplying  by  g/U  and  noting  T 

m 

=  t  following  nondimensional  expression  results  which  may  be  compared 

to  Bretschneider ’ s  result: 


U  -  0.9 


(34) 


This  comparison  is  made  in  Figure  5,  and  it  demonstrates  a  reasonable  fit  to 
the  field  data.  The  advantage  of  Equation  33  is  that  a  value  of  windspeed  is 
not  a  factor  in  determining  the  fully  developed  cutoff  frequency  in  shallow 
water.  In  otjier  words,  it  is  assumed  that  once  a  fully  developed  condition 
exists  in ^shallow  water,  the  peak  frequency  remains  fixed  regardless  of  any 
increase  in  windspeed. 
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Figure  5.  Relationship  between  dimensionless  period 
and  depth  (dashed  line  is  Equation  34;  solid  line  is 
Bretschneider ' s  (1958)  empirical  fit)  (Vincent  and 

Hughes  1985) 


Wave  height  limit 

30.  Bretschneider  (1958)  also  presented  nondimensional  significant  wave 
height  versus  nondimensional  depth  for  the  shallow-water  fully  developed  data. 
Vincent  and  Hughes  (1985)  were  able  to  show  that  Vincent's  estimate  for  depth- 
limited  significant  wave  height  (Equation  20)  can  be  used  to  derive  an  expres¬ 
sion  for  the  fully  developed  case  closely  matching  Bretschneider ' s  empirical 
formulation.  They  expressed  in  Equation  24  as 


k  = 
m 


L 

m 


271 


(gh)^^^T 


m 


0.9 

h 


where  is  replaced  by  its  shallow~water  expression  from  linear  theory  and 

is  found  from  Equation  33  for  the  fully  developed  case.  Using  this  expres¬ 
sion  for  ,  the  value  of  a  for  the  fully  developed  case  is  determined 
from  Equation  22  as 

a  =  0.0074U/(gh)^^^  (35) 


where  the  0.49  exponent  of  Equation  22  is  taken  as  1/2.  Substituting  Equa¬ 
tion  35  for  a  and  Equation  33  for  into  Equation  20  yields 


«£  = 


0.210U^^^h^^^ 


(36) 


where  H.  is  the  significant  wave  height  for  the  fully  developed  case.  Non- 

^  2 
dimensionalizing  by  g/U  gives 


gH, 


U‘ 


=  0.210 


(37) 


Equation  37  is  compared  to  Bretschneider' s  result  in  Figure  6, 

31.  The  favorable  fit  of  the  equations  for  fully  developed  wave  period 
and  wave  height  to  the  shallow-water  field  data  of  Bretschneider  lend  support 
to  the  TMA  shallow-water  self-similar  spectral  form  and  to  the  parameteriza- 
tions  made  for  a  .  However,  there  are  still  many  unanswered  questions  in¬ 
volving  the  underlying  physical  processes  which  result  in  the  TMA  spectrum. 

It  must  be  noted  that  these  shallow-water  limitations  are  for  full  development 
on  a  flat  bottom,  but  they  do  provide  a  useful  upper  limit  for  some  engineer¬ 
ing  purposes . 
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Figure  6.  Relationship  between  dimensionless  wave 
height  and  depth  (dashed  line  is  Equation  37; 
solid  line  is  Bretschneider* s  (1958)  empirical  fit 
(Vincent  and  Hughes  1985) 


PART  III:  SURF  ZONE  SPECTRA 


32.  In  the  formulation  and  parameterization  of  the  TMA  spectral  form 
by  Bouws  et  al.  (1985a)  the  data  examined  were  restricted  to  those  spectra  re¬ 
corded  outside  of  the  surf  zone.  They  were  uncertain  of  the  appropriateness 
of  using  linear  wave  theory  in  the  surf  zone  to  convert  from  wave  number  space 
to  frequency  space  via  Kitaigorodskii  et  al.'s  (1975)  <l>  function.  Also, 

outside  the  surf  zone  the  balance  of  forces  determining  spectral  evolution  are 
comprised  of  wind,  breaking  (largely  at  high  frequencies),  dissipation  due  to 
a  bottom  boundary  layer,  and  transfers  of  energy  in  the  spectrum.  Once  the 
surf  zone  is  reached,  major  breaking  occurs  near  the  peak  of  the  spectrum  thus 
introducing  a  separate  concern.  Vincent  examined  the  fit  of  the  TMA  spectral 
form  to  spectra  collected  well  within  the  surf  zone  during  the  ARSLOE  experi¬ 
ment.  (The  surf  zone  spectra  were  not  a  part  of  the  ARSLOE  portion  of  the 
data  set  used  to  determine  the  TMA  parameters).* 


Fit  of  TMA  to  Collected  Spectra 


33.  The  TMA  form  was  fitted  to  the  spectra  only  over  the  region  up  to 

twice  the  spectral  peak  frequency  2f  .  Beyond  this  frequency  forced  har- 

m 

monies  begin  to  dominate  the  spectrum,  while  the  TMA  form  is  meant  to  represent 

only  free  wave  energy.  Figure  7  is  an  example  of  the  best  fit  obtained  by 

Vincent  to  a  surf  zone  spectrum  recorded  in  1.7  m  of  water.  A  reasonable  fit 

is  mace  between  the  range  of  0.5f  to  2. Of  ,  which  is  the  principal  energy 

m  m 

containing  region.  At  higher  frequencies  the  TMA  form  underestimates  the 
energy  content.  Vincent  (1985b)  noted  that  the  energy  contained  at  these 
higher  frequencies  increases  as  the  depth  decreases.  This  increased  energy 
represents  the  increasing  nonlinearity  of  the  long  period  wave  components  in 
shallow  water.  Vincent  suggested  a  modifying  factor  for  the  TMA  spectral  form, 
but  he  lacked  the  necessary  data  to  provide  a  generally  applicable  modification. 
However,  as  a  first  approximation,  the  unmodified  TMA  spectrum  can  be  used  to 
represent  surf  zone  spectra  if  the  spectral  parameters  can  be  estimated.  The 
resulting  estimated  significant  wave  height  found  using  this  approach  was 
within  10  percent  of  the  actual  values  for  the  cases  that  Vincent  examined. 

*  Personal  correspondence  between  Vincent  and  Hughes,  January  1984. 
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Figure  7.  Fit  of  TMA  spectrum  to  surf  zone 
irregular  wave  spectrum 

l.ote  that  the  logarithmic  scale  for  energy  density  used  in  Figure  7  makes  it 
appear  that  the  difference  in  total  energy  between  the  TMA  fit  and  the  actual 
spectrum  is  greater  than  it  truly  is.) 

34.  An  unexpected  result  of  this  work  was  the  reasonable  fit  of  the 
TMA  spectral  form  to  broad  swell  spectra  in  the  surf  zone  obtained  from  both 
the  field  and  the  laboratory.  This  is  surprising  for  two  reasons.  First, 
the  broad  swell  is  more  indicative  of  a  decaying  than  a  growing  sea  state  for 
which  the  self-similar  spectrum  was  proposed.  Second,  there  is  no  reason  to 
expect  that  mechanically  generated  irregular  waves  in  the  laboratory  should 
follow  the  TMA  form  for  wind  waves  since  the  primary  forcing  function  (the 
wind)  is  absent.  This  latter  observation  was  found  to  be  true  for  both  char¬ 
acteristic  broad  swell  spectra  outside  the  surf  zone  as  well  as  for  typical 
wind  sea  spectra  simulated  in  the  laboratory.  Figure  8  illustrates  the  fit¬ 
ting  of  the  TMA  spectral  form  to  laboratory  data  shoaled  over  a  1:30  slope. 

Hypothesis  for  Energy  Saturation 

Local  saturation  zone 

35.  Vincent  (1985b),  suggested  a  hypothesis  for  energy  saturation  of 
irregular  waves  during  shoaling  by  proposing  a  separation  of  the  saturation 
process  during  shoaling  into  two  distinct  types:  the  local  and  universal 

20 
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saturation  zones.  As  waves  progress  into  finite  depth  waters,  they  first 
enter  a  local  saturation  zone  in  the  intermediate  and  shallow-water  depths 
where  nonlinear  wave-wave  interactions  shift  energy  away  from  the  principal 
energy  containing  frequencies  to  higher  frequencies  where  it  is  dissipated  by 
deepwater  breaking  mechanisms.  This  is  the  region  for  which  all  the  previous 
TMA  formulations  and  parameterizations  for  wind  waves  are  valid.  It  appears 
that  broad  swell  spectra  also  can  be  fitted  by  the  TMA  form  in  this  region, 
but  there  is  no  external  parameterization  for  a  in  this  case. 

Universal  saturation  zone 

36.  Eventually,  if  the  water  depth  becomes  sufficiently  shallow,  wave 
breaking  begins  to  occur  at  the  peak  of  the  spectrum,  and  dissipation  dominates 
the  spectrum.  Vincent  (1985b)  referred  to  this  region  as  the  universal  satura¬ 
tion  zone.  While  there  still  may  be  some  transfer  of  energy  between  frequen¬ 
cies  in  this  zone,  the  wave  train  is  basically  a  series  of  nondispersive  waves, 
each  breaking  in  sequence.  The  fact  that  the  TMA  spectral  form  could  be  fitted 
to  data  gathered  within  the  surf  zone  indicated  to  Vincent  that  the  spectral 
parameters  could  be  obtained  from  the  steepness  equations  (Equations  26  and  28) 
if  a  saturated  value  for  the  energy  could  be  obtained.  Working  on  the  hypothe¬ 
sis  that  the  largest  wave  at  a  given  depth  in  the  universal  saturation  zone  is 
controlled  by  the  depth,  Vincent  (1985b)  proposed 


H  =  Bh 
mo 

where 

h  = 

water  depth 

B  = 

R 

max 

R  R 
so  ms 

and 

R  =  H  /h 

max  max 

for  monochromatic  waves 

R  =  H  /H 

so  s  mo 

where 


(38) 


(39) 


H  =  average  of  the  highest  1/3  waves 


Values  for  can  be  obtained  from  Figure  9,  which  is  reproduced  from  the 

SPM  (1984),  and  values  for  can  be  found  from  the  "average"  line  in  Fig¬ 

ure  10,  which  is  reproduced  from  Thompson  and  Vincent  (1985).  For  slopes 
typical  of  most  beaches,  the  coefficient  B  ranges  between  0.55  and  0.65. 

These  values  are  consistent  with  the  value  of  B  =  0.6  suggested  by  Thornton 
and  Guza  (1982),  and  unless  the  area  of  concern  varies  greatly  from  a  "typical" 
beach,  the  value  of  0.6  for  B  is  probably  sufficient  for  engineering  pur¬ 
poses.  The  significant  steepness  is  found  from  Equation  27  as 


e  = 


_  (E)^^^  "mo 


m 


H 

mo 

4L 

m 


Bh 


4(gh)^^^T 


m 


(40) 


when  Equation  19  is  used  for  and  the  linear  theory  shallow-water  wave¬ 

length  is  substituted.  The  parameter  a  can  now  be  found  from  Equation  26 


where 


(41) 


(42) 


37.  The  ability  of  Equation  41  to' predict  a  in  very  shallow  water 
shown  in  Figure  11  (from  Vincent  1984).  The  figure  shows  reasonable  agree¬ 
ment  for  both  laboratory  and  field  data  in  the  universal  saturation  zone. 

Note  that  in  the  universal  saturation  zone  H  varies  linearly  with  depth 

mo 

(Equation  38) ,  while  in  the  local  saturation  zone  varies  with  the  square 

root  of  the  depth  (Equation  20).  Also  note  that  the  equation  for  in  the 

local  saturation  zone  (Equation  20)  reduces  to  the  expression  for  in  the 

I 

universal  saturation  zone  (Equation  38)  when  Equation  41  is  substituted  for 
a  .  This  is  not  surprising  since  both  equations  are  linked  to  the  equation  for 
a  in  terms  of  e  . 

38.  It  is  somewhat  difficult  to  define  a  demarcation  line  between  the 
local  saturation  zone  and  the  universal  saturation  zone.  Ideally  there  should 
be  some  limiting  value  of  significant  steepness  beyond  which  wave  breaking  at 
the  spectral  peak  is  assured. 
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Figure  9.  Dimensionless  breaker  height  versus  relative  depth  (from  SPM  1984) 


PART  IV:  DISCUSSION 


39.  Thus  far  the  focus  of  this  report  has  been  on  past  developments  in 
the  specification  of  self-similar  spectral  forms  and  their  role  in  the  develop¬ 
ment  of  the  new  TMA  spectral  form.  By  adding  more  and  more  terms  to  Phillips' 
(1958)  representation  of  the  limiting  energy  density  in  the  equilibrium  range, 
researchers  progressively  specified  a  self-similar  spectral  shape  for  deep 
water  (Pierson-Moskowitz) ,  a  fetch-limited  deepwater  spectrum  (JONSWAP) ,  a 
shallow-water  equilibrium  range  (Kitaigorodskii  et  al,),  and  finally  a  self¬ 
similar  equilibrium  spectral  form  valid  at  all  water  depths  (TMA) . 

40.  The  self-similar  parametric  TMA  representation  was  shown  to  have 
the  ability  to  represent  wind  waves  and  broad  swell,  and  it  appears  that  both 
the  deepwater  and  shallow-water  limits  merge  with  previously  obtained  empirical 
results.  It  was  seen  also  that  the  TMA  form  even  fit  very  shallow  surf  zone 
spectra  up  to  twice  the  peak  frequency. 


Physical  Implications  of  the  TMA  Spectral  Form 


41.  The  fact  that  the  self-similar  TMA  form  can  describe  wind  sea  spec¬ 
tra  in  deep,  intermediate,  or  shallow  water  reinforces  the  conclusion  of 
Kitaigorodskii  et  al.  (1975)  that  the  basic  scaling  lies  in  k-space  rather 
than  in  frequency  space.  Also,  the  linear  theory  dispersion  relation  seems 
adequate  for  converting  from  wave  number  space  to  frequency  space  via  the  <l> 
factor. 

-3 

42.  Both  the  JONSWAP  and  the  TMA  spectral  forms  extended  the  k  scal¬ 
ing  for  the  equilibrium  range  across  the  entire  spectrum.  There  is  no  physical 
explanation  as  to  why  this  scaling  should  be  valid,  but  years  of  experience 
using  the  JONSWAP  form,  coupled  with  the  relative  good  fit  of  the  TMA  form  to 
shallow-water  spectra,  indicates  that  the  scaling  does  hold  outside  the  equi¬ 
librium  range. 

43.  Self-similarity  in  spectral  form  implies  that  there  is  a  relative 
balance  between  the  wind  energy  input,  energy  transfers  within  the  spectrum, 
and  energy  dissipation.  This  balance  is  maintained  and  limited  to  produce  a 
consistent  spectral  shape.  No  analyses  have  been  made  into  the  very  difficult 
problem  of  the  spectral  balance  between  the  various  source/sink  terms  such  as 
nonlinear  wave-wave  interactions,  wave  dissipation,  atmospheric  input,  bottom 
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friction,  percolation,  etc.,  for  the  TMA  data  set.  However,  it  is  tentatively 
concluded  that  the  bottom  friction  effect  in  the  dissipation  of  wave  energy  is 
substantially  less  than  originally  thought.  This  conclusion  is  based  on  the 
observation  that  the  parametric  relationships  provide  equally  good  spectral 
approximations  over  a  wide  variety  of  depths  and  bottom  materials.  Thus  it  is 
unlikely  that  bottom  dissipation  would  be  a  dominant  influence  on  spectral 
shape.  With  the  establishment  of  the  TMA  shape,  study  into  the  various  source/ 
sink  terms  now  can  be  restricted  to  combinations  which  will  produce  a  self¬ 
similar  shape.  Bottom  friction  for  typical  beaches  is  thought  to  be  more  im¬ 
portant  for  swell  than  for  wind  seas. 

44.  Kitaigorodskii  et  al.’s  (1975)  equilibrium  range  was  intended  for 
the  case  of  a  horizontal  bottom,  but  the  parameterization  and  application  of 
the  TMA  spectrum  to  data  from  gently  sloping  bottoms  up  to  slopes  of  1:100 
indicate  that  the  self-similar  equilibrium  form  evolves  rapidly  enough  to  ade¬ 
quately  describe  the  spectral  transformation  over  gentle  slopes.  Steep  sloped 
beaches  might  not  permit  the  equilibrium  condition  to  be  reached;  hence,  the 
TMA  spectrum  could  well  represent  an  unattainable  limit  on  a  steep  slope. 

45.  Vincent  (1985b)  sums  up  the  TMA  spectral  form  as  "...a  quantitative 
description  of  the  wave  transformation  process  with  a  useful  parameterization 
scheme."  He  adds  that  while  many  of  the  principal  physical  processes  occurring 
in  wave  shoaling  have  been  identified,  a  deterministic  description  is  still  in 
the  future. 


Jtions  and  Limitations  of  the  TMA  Spectral  Form 


46.  The  primary  underlying  assumption  invoked  when  applying  the  TMA 
spectrum  and  the  prognostic  equations  is  that  the  wind  sea  is  at  a  steady 
state  or  equilibrium  condition.  This  means  that  the  wind  has  been  steady  long 
enough  (as  yet  undefined)  for  the  waves  to  reach  equilibrium  and  that  the 
bottom  topography  is  a  gentle  slope  with  smoothly  varying  features  and  without 
complexities  which  might  cause  rapid  alteration  of  the  wave  train.  Until 


further  data  are  examined,  a  slope  steepness  guideline  of  1:100  is  suggested 
which  arises  from  the  limit  of  the  TMA  data  set.  While  steeper  slopes  might 
be  able  to  maintain  the  TMA  equilibrium  condition,  the  worst  that  can  happen 


by  using  the  TMA  approximation  on  a  steep  slope  is  a  conservative  estimate  of 
the  maximum  possible  significant  wave  height  for  the  given  conditions.  Note 
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that  laboratory  data  on  a  1:30  slope  were  adequately  described  using  the  TMA 
representation. 

47 .  The  TMA  spectrum  as  it  is  presently  parameterized  cannot  be  used 
in  fetch-  or  duration-limited  shallow-water  wave  growth  situations  since  it 
is  a  final  steady  state  form.  Fetch-  or  duration-limited  design  in  shallow 
water  needs  to  be  conducted  using  numerical  models  (for  complex  regions)  or 
wave  forecasting  curves  such  as  provided  in  ETL  1110-2-305  (Vincent  and  Lock¬ 
hart  1984).  Future  developments  will  probably  provide  a  TMA  fetch-  and 
duration-limited  form  paralleling  the  deepwater  JONSWAP  equations. 

48.  The  TMA  spectral  form  is  for  a  single-peaked  spectrum  containing 
only  minor  low-frequency  energy.  In  most  instances  large  wave  conditions  will 
result  from  either  a  local  storm  or  swell  resulting  from  a  nearby  storm.  These 
spectra  will  probably  propagate  toward  shore  as  a  single-peaked  spectrum  de- 
scribable  by  the  TMA  shape.  The  TMA  formulation  does  not  handle  multipeaked 
spectra  which,  according  to  Thompson  (1980),  constitute  over  60  percent  of  the 
coastal  wave  conditions  in  the  United  States.  However,  most  of  the  multi¬ 
peaked  cases  are  of  low  energy  and  of  minor  importance  for  engineering  design. 
The  TMA  shape  may  be  used  to  describe  the  individual  spectral  peaks. 

49.  The  shallow-water  limits  to  wave  growth  derived  as  a  cutoff  fre¬ 
quency  and  a  limiting  significant  wave  height  fit  the  field  data  and  provide 
useful  engineering  guidance.  The  limits  are  strictly  for  full  development 
over  a  horizontal  bottom  or  a  gently  sloping  bottom. 

50.  The  TMA  prognostic  equations  require  specifications  of  either 

(a)  U  ,  h  ,  and  f  ,  or  (b)  e  =  f(H  ,  h  ,  and  f  )  .  Vincent  (1984) 

m  mo  m 

suggests  that  when  wave  conditions  are  steady  and  the  combination  of  beach 

slope  and  propagation  distance  are  sufficiently  small,  then  f  may  be  pro- 

m 

jected  from  deep  water  into  shallow  water.  In  other  words,  the  peak  frequency 
doesn't  shift  during  shoaling.  Field  data  from  two  1980  storms  at  Duck,  North 
Carolina,  indicated  that  f^  did  not  change  more  than  10  percent  over  a  36  km 
distance  in  which  the  depth  changed  from  36  to  2  m.  In  cases  where  wind 
fields  and  bottom  conditions  are  highly  inhomogeneous  and  f^  is  expected  to 
shift,  a  fully  time-dependent  spectral  model  is  needed  to  determine  f^  in 
shallow  water. 
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PART  V:  ENGINEERING  APPLICATIONS 


51.  The  preceding  formulations  for  a  self-similar  equilibrium  spectrum 
for  all  water  depths  can  provide  useful  engineering  design  guidance  because 
the  equilibrium  spectrum  puts  constraints  on  the  expected  energy  levels  and 
spectral  shapes  under  given  conditions.  While  it  is  impossible  to  list  every 
conceivable  application  of  the  TMA  spectrum,  the  following  short  examples  il¬ 
lustrate  several  engineering  uses. 

Numerical  Models 

52.  Some  of  the  current  shallow-water  wave  growth  and  propagation 
numerical  models  can  simulate  wave  propagation  over  complex  bathymetry  during 
unsteady  wind  conditions.  These  models  typically  contain  many  source/sink 
terms,  some  of  which  can  be  quite  site  specific.  In  some  applications  it  is 
possible  for  the  models  to  give  exaggerated  results  because  there  is  no  effec¬ 
tive  upper  limit  to  the  wave  growth.  This  can  be  corrected  by  constraining 
the  numerical  simulation  to  an  upper  limit  of  growth  on  the  high  frequency 
side  of  the  spectral  peak  as  defined  by  the  TMA  spectrum. 

53.  Oceanographers  at  CERC  illustrated  the  usefulness  of  the  TMA  spec¬ 
tral  form  as  an  upper  limit  when  they  modified  an  existing  numerical  wave 
growth  spectral  model  and  ran  several  simulations,  with  and  without  the  TMA 
prescribed  limit.  Generally,  it  was  found  that  this  particular  numerical 
model  yielded  unrealistically  high  results  for  high  windspeeds  and  for  beach 
slopes  of  1:200.  In  other  words,  sufficient  energy  dissipation  through  bottom 
friction  and  other  sink  terms  did  not  occur.  When  the  slope  was  very  mild 
(1:2,000),  bottom  friction  attenuated  the  waves  adequately.  Incorporating  the 
TMA  limiting  spectral  form  into  the  numerical  model  has  made  it  a  better  en¬ 
gineering  tool;  however,  this  does  not  relieve  the  engineer  from  responsible 
calibration  and  verification  of  the  model. 

Design  Example  1 


Given 

54.  For  a  smoothly  varying  bottom  slope,  construct  a  set  of  design 

curves  showing  the  depth-limited  significant  wave  height  H  as  a  function 

mo 


of  peak  spectral  period  and  windspeed  U  .  The  given  water  depth  is  5  m 

(16  ft). 

Solution 

55.  It  is  first  necessary  to  choose  the  proper  equation  for  H 

mo 

(either  Equation  29  or  30).  This  is  done  by  determining  the  wave  frequency 
where  1  »  i -e. , 


=  27tf(h/g)^/^  =  1 


or 


f’  =  (9.81/5)^^^/2n  =  0.22  Hz 


which  is  T*  a«  4.5  sec.  When  the  principal  energy  containing  frequencies 

have  periods  greater  than  4.5  sec.  Equation  30  can  be  used.  We  will  assume 

that  f  =  1.5  f  or  T  =  1.5  T'  =  6.75  sec.  Any  calculations  made  for 
mm 

peak  periods  less  than  6.75  sec  should  be  done  using  Equation  29  with 
being  determined  using  the  intermediate  depth  linear  dispersion  relation. 
The  calculations  are  done  as  follows: 


For  U  =  10  m/sec  (32.8  ft/sec)  and  T  =  10  sec,  first  find 


m 


=  (gh)^^^T^  =  [(9.81)(5)]^^^(10)  =  70  m  (230  ft) 


Next  find  K  from  Equation  24. 


K  = 


2 

U  2n 
*  m 


QO)  (2n)  ^  0  9J5 

(9.81)(70) 


56.  Using  Equation  22, 


a  =  0.0078  =  0.0078  (0.915)°’^®  =  0.00747 


and  the  H  associated  with  this  T  and  U  is  by  Equation  30 
mo  m  ^ 

H  =  I  (08h)^^^T  =  I  r(0.00747)(9.8l)(5)V^^  •  (10)  =  1.93  m  (6.3  ft) 

mo  71  m  IT  L  J 

To  be  slightly  conservative  (10  percent),  use  the  1.1  factor  as  given  in 
Vincent's  original  derivation  (Equation  20)  to  finally  arrive  at 
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H  =  1.1  (1.93)  =  2.1  m  (6.9  ft) 
mo 

A  simple  programming  of  these  few  steps  supplied  the  repetitive  calculations 
necessary  to  produce  the  set  of  curves  given  in  Figure  12  for  the  stated  ini¬ 
tial  conditions.  The  region  of  breaking  is  given  by  Equation  38  as 

H  =Bh=0.6(5m)=3.0m 
mo 


Design  Example  2 

Given 

57.  To  determine  crest  elevation  for  some  structures  where  limited  over 
topping  is  allowable,  engineers  should  use  the  following  wave  statistics:  H 

s 

(average  of  the  highest  1/3  waves),  (average  of  the  highest  10  percent 

of  all  waves),  and  (average  of  the  highest  1  percent  of  all  waves). 

They  should  then  calculate  these  values  for  a  windspeed  of  10  m/sec  (32.8  ft/ 
sec)  and  a  peak  period  of  10  sec  in  a  water  depth  of  5  m  (16.4  ft).  (Note 
that  these  are  the  same  values  used  in  the  sample  calculation  of  Example  1.) 
Solution 

58.  From  Example  1  it  is  seen  that  the  given  values  result  in  = 

2.1  m  (6.9  ft).  First  determine  if  this  is  a  breaking  condition  as  defined 

by  H  =  Bh  when  B  =  0 . 6  ,  i . e . , 
mo 

H  ,,  ,  .  =  0.6  (5  m)  =  3.0  m 

mo (break) 


Since  =  2.1  m  is  nonbreaking,  we  can  determine  from  Figure  10  using 


The  values  for  and 

10% 

the  SPM  (1984). 


are  now  found  using  the  expressions  given  in 


=  1.27  H 
10%  s 


H 


1% 


1.67  H 

s 


1.27(2.48)  =  3.15  m  (round  off  to  3.2  m)  (10.3  ft) 
1.67(2.48)  =  4.14  m  (round  off  to  4.2  m)  (13.6  ft) 


Note  that  the  above  calculations  imply  that  the  Rayleigh  distribution  of  wave 
heights  is  valid  under  the  given  conditions,  while  the  highest  possible  wave 
at  this  depth  by  monochromatic  theory  is  about  0.76  (5  m)  =  3.8  m.  Thus  the 
tail  of  the  wave  height  distribution  will  be  truncated  by  breaking. 


Design  Example  3 


Given 

59.  Determine  the  conditions  under  which  a  project  built  in  7  m  (23  ft) 
of  water  would  be  subjected  to  breaking  wave  conditions  (i.e.,  surf  zone  con¬ 
ditions)  caused  by  irregular  waves  with  a  peak  period  of  11  sec. 

Solution 

60.  The  for  the  breaking  limit  (universal  saturation)  is  found 

from  Equation  38  as 


H  =  Bh  =  0.6  (7  m)  =  4.2  m  (13.8  ft) 
mo 


The  spectral  parameter  a  can  be  found  from  Equation  41  with  tUj^  determined 
from  Equation  42,  i.e.. 


hm 


a 


\l/2 


(i)  =  TT  (9^) 

^Bu.j^/2^^  =  [(0.6)(0.4825)/2] 


0.4825 

^  =  0.021 


61.  Next,  use  Equation  22  to  find  what  value  of  K  will  produce  the 
following  value  of  a  at  the  given  depth  and  peak  period: 


0.49 


a  =  0.0078k 


or 


K  =  (a/0. 0078)^^°*^^  =  (0. 021/0. =  7.55 


From  Equation  24  with 


,  2n  2n 
k  =  —  = 
m  L 

m 
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(gh)^''^T^  [(9.81)(7)]^/^(11) 


-  0.0689  rad/m 


we  get 


U  =  (xg/k^y^^  =  [(7.55)(9.81)/0.689]^^^  =  32.8  m/sec  (107.6  ft/sec) 


This  windspeed  is  the  sustained  wind  necessary  to  develop  the  equilibrium 
breaking  condition  at  the  site  for  the  given  peak  period.  It  corresponds 
to  speeds  on  the  order  of  73  mph  which  is  a  full  hurricane  force  wind. 


Design  Example  4 


‘*r-v 


Given 

62.  A  wind  sea  has  a  significant  wave  height  of  2  m  and  a  peak 

spectral  period  of  8  sec  in  a  water  depth  of  10  m.  Determine  the  corre¬ 
sponding  in  4-m  water  depth,  assuming  no  refraction  or  diffraction  and 

that  the  peak  period  remains  constant. 

Solution 

63.  From  Equation  29, 

Let 

in  10  m  in  4  m 


H  =  H 
1  mo 


Oj  =  a 


L  =  L 
1  m 


H_  =  H 
2  mo 


a^  =  a 


^9  =  K 
2  m 


Taking  the  ratio  of  these  two  relationsnips  yields 
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1/2 


The  ratio  using  Equation  22  as 


Substituting  this  relationship  into  the  previous  equation  gives 


which  is  good  under  the  restrictive  assumptions  of  no  refraction  and  diffrac¬ 
tion  and  that  the  wind  is  equal  and  constant  at  both  sites. 

65.  Using  linear  wave  theory  to  find  the  wavelength  associated  with  the 
peak  spectral  period  (assumed  to  remain  constant  at  8  sec): 


Lj  =  70.9  m  at  h  =  10  m 


L2  =  48.0  m  at  h  =  A  m 


H„  \48.0/ 


=  1.34 


so  that 


T =1.5m  =  H  at  a  water  depth  h  =  4  m 
1.34  mo 


If  the  significant  wave  height  had  been  shoaled  using  linear  wave  theory, 


H,/H-  would  have  been  0.7  and  H  at  h  =  4  m  would  have  been  2.9  m, 
12  mo 


PART  VI :  SUMMARY 


66.  Recent  work  by  Bouws  et  al,  has  led  to  a  self-similar  equilibrium 
spectral  form  for  all  water  depths.  This  form,  called  the  TMA  spectrum,  arises 
from  the  substitution  of  Kitaigorodskii  et  al.'s  finite  depth  equilibrium  range 
for  the  deepwater  equilibrium  range  in  the  JONSWAP  energy  density  equation. 
Examination  of  over  2,800  field  spectra  in  varying  water  depths  and  under  di¬ 
verse  conditions  resulted  in  a  parameterization  for  the  spectral  variables  in 
terms  of  depth,  peak  frequency,  and  windspeed.  The  parameters  Of  and  y 

were  also  expressed  in  terms  of  significant  wave  steepness.  A  simple  expres¬ 
sion  for  the  energy-based  significant  wave  height  H  was  derived  from  the 

mo 

TMA  spectral  form  making  it  possible  to  predict  the  depth-limited  equilibrium 

H  in  either  intermediate  or  shallow  water, 
mo 

67.  The  TMA  spectral  representation  recovers  the  deepwater  fully  devel¬ 
oped  results  when  the  deepwater  limits  are  applied,  and  it  follows  the  fully 
developed  shallow-water  empirical  results  when  the  suggested  cutoff  frequency 
relationship  for  shallow  water  is  applied. 

68.  Within  the  surf  zone  the  TMA  spectrum  can  be  fit  rather  well  up  to 
twice  the  peak  frequency.  Beyond  that  point  highly  nonlinear  processes  cause 
a  deviation  from  the  TMA  form.  A  distinction  was  made  between  local  and  uni¬ 
versal  saturation,  and  a  method  of  estimating  the  spectral  parameters  within 
the  surf  zone  was  presented. 

69.  The  success  of  the  TMA  spectrum  as  a  simple  spectral  model  rein¬ 
forces  Kitaigorodskii  et  al.'s  (1975)  conclusion  that  the  basic  self-similarity 
scaling  is  in  k-space.  Many  of  the  principal  physical  processes  which  bring 
about  the  TMA  form  have  been  identified  but  not  quantified.  Applications  and 
examples  were  presented  to  illustrate  some  uses  of  the  TMA  results  under  the 
stated  assumptions  and  limitations. 
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Parameter 

Relative  depth 

Total  spectral  energy 

Energy  density  for  JONSWAP  spectrum 

Equilibrium  energy  density  in  frequency  space 

Finite  depth  equilibrium  energy  density  in  frequency  space 

Energy  density  for  Pierson-Moskowitz  spectrum 

Energy  density  for  TMA  spectrum 

Frequency 

Frequency 

Highest  freqeuncy  where  shallow-water  dispersion  relation  holds 
Cut  off  frequency 

Equilibrium  energy  density  in  wave  number  space 
Spectral  peak  frequency 
Gravitational  constant 
Water  depth 

H  for  fully  developed  wind  sea  over  a  flat  bottom 
mo 

Maximum  monochromatic  wave  height 
Zero-moment  wave  height 
Average  height  of  the  1/3  highest  waves 
Average  of  highest  1  percent  of  waves 
Average  of  highest  10  percent  of  waves 
Wave  number 

Wave  number  for  waves  at  peak  frequency 
Wavelength 

Wavelength  associated  with  frequency  f^  by  linear  wave  theory 
Deepwater  equivalent  of 

Relative  maximum  wave  height  for  monochromatic  waves 

Ratio  of  H  to  H 
max  s 

Ratio  of  H  to  H 

s  mo 

Period 

Spectral  peak  period 
10-m-high  windspeed 
Fetch  distance 
Spectral  parameter 
A  constant 


6 


Spectral  parameter 
Significant  wave  steepness 
Dimensionless  wave  number 

3.14... 

Spectral  parameter 
Spectral  parameter 
Spectral  parameter 
Finite  depth  factor 
Circular  frequency 
Depth  dependent  frequency 
(ii.  evaluated  at  f_ 


